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Within a path integral formalism for non-Gaussian price fluctuations we set up a 
simple stochastic calculus and derive a natural martingale for option pricing from the 
wealth balance of options, stocks, and bonds. The resulting formula is evaluated for 
truncated Levy distributions. 

^ " I. INTRODUCTION 

O ; 

Financial markets exhibit non-Gaussian price fluctuations. For this reason they are incomplete markets 
and mathematicians are unable determine a unique riskless martingale distribution to calculate fair option 
, prices. 

^ ■ As pointed out in Pareto's 19th century work |jj and reemphasized by Mandelbrot in the 1960s, the 
logarithms of assets prices in financial markets do not fluctuate with Gaussian distributions, but possess 
0^ ' much larger tails which may be approximated by various other distributions such as truncated Levy 
distributions , Meixner distributions HfjJ , generalized hyperbolic distributions or simplifies versions 
thereof This has the unpleasant consequence that the associated stochastic differential equations 

cannot be treated with the popular Ito calculus. The mathematicians call such markets incomplete, 
implying that there are different choices of riskless martingale distributions with which one can calculate 
I ■ option prices. Many of these have been discussed in the literature, and mathematicians have invented 
various sophisticated criteria under which one would be preferable over the others for calculating financial 
CN| ■ risks. Davis, for instance, has introduced a so-called utility function p3[ which is supposed to select 
| optimal martingales for different purposes. 

In this paper we want to point out that in a path integral formulation of the problem, a straightforward 
extension of the old chain of arguments which led Black-Scholes to their famous formula produces a specific 
simple martingale which is different from presently popular versions based on Esscher transforms. Before 
we come to this we briefly review the derivation of a stochastic calculus for non-Gaussian processes |34| 
from path integrals which extends Ito's calculus in the simplest way. 

c . 

O . II- GAUSSIAN APPROXIMATION TO FLUCTUATION PROPERTIES OF STOCK PRICES 

O 

Let S(t) denote the price of some stock. On the average, stocks grow exponentially, and a first 
approximation considers a stock price as an exponential S(t) = e xs ^' of a fluctuating variable x(t) which 
obeys a stochastic differential equation 



xs(t) =r xs + a7?(i), (1) 
where rj{t) is a white noise of unit strength with the correlation functions 

fa(t)>=0, ( V (t)r 1 (t')) = S(t-t'). (2) 
The stock price S(t) itself satisfies the stochastic differential equation for exponential growth 

S(t) 



S(t) 



rs + (TT)(t), (3) 



where The standard deviation a is a measure for the volatility of the stock price, which defined by the 
expectation value 
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'S(jt) 




S(t) 



(4) 

The rate constants in Eqs. ([!]) and (||) differ from one another due the stochastic nature of xs(t) and 
S(t). According to Ito's rule, we may expand 



dx s fy, 



1 d 2 x s ■ 



dS 
S(t) 



2 dS 2 
S(t) 



S(t) 



S\t)dt 



dt - 



(5) 



and replace the last term by its expectation value 

• / \ S(t) 
±s(t) - 




dt 



S(t) 2 (6) 
Inserting Eq. (||), we obtain the well-known Ito relation 

r xs = r s - \f 2 . (7) 

In praxis, this relation implies that if we fit a straight line through a plot of the logarithms of stock prices, 
the forward extrapolation of the average stock price is given by 



S(t) = S(0)e 



rgt 



S(0)e 



(r xs +<T*/2)t 



(8) 



III. NON-GAUSSIAN DISTRIBUTIONS 



The description of the fluctuations of the logarithms of the stock prices around the linear trend by a 
Gaussian distribution is only a rough approximation to the real stock prices. Given a certain time scale, 
for instance days, they follow some distribution D(x) which is the Fourier transform of some function 
D(p), which we shall write as an exponential e~ H ^: 



J 2-7T J 2n 

In general, H (p) will have a power series expansion 



H(p) = ic x p+ 7;c 2 p 2 - «^j c 3P 3 " ^ C4p4 + i ^ C5p5 



(9) 



(10) 



The coefficients c n in the expansion (|T^) are the cumulants of the distribution D(x), from which they 
can be obtained by the connected expectation values 

c\ = (x) c = (x) = dxx D(x), C2 = (x 2 ) c = (x 2 ) — (x) 2 = / dx (x — (x)) 2 = a 2 D(x), . . . 



By analogy with mechanics, we call H (p) the Hamiltonian of the fluctuations. A Gaussian Hamiltonian 
H{p) = a 2 p 2 /2 coincides with the mechanical energy of a free particle of mass 1/er 2 . Non-Gaussian 
Hamiltonians are standard in elementary particle physics. A relativistic particle, for example, has a 
Hamiltonian H(p) — \/p 2 + m 2 , and there is no problem in defining and solving the associated path 
integral p5| |. 

We also introduce a modified Hamiltonian H(p), whose Fourier transform has zero average. Its expan- 
sion ( [Toj ) has the linear term ic\p subtracted. We also define a modified Hamiltonian H r {p) by adding to 
H (p) a linear term irp, i.e., 



H(p) = H{p) - H'(0)p, H r {p) = H{p) + irp = H{p) - H'(0)p - 



irp. 



(12) 
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A. Path Integral for Non-Gaussian Fluctuations 



It is easy to calculate the properties of a process (Q) whose fluctuations are distributed according to 
a general non-Gaussian distribution. If we assume the rate r xs to coincide with the linear coefficient c\ 
in the expansion (|l0|), such that we can abbreviate -ffr xs (p) by H(p), the stochastic differential equation 
reads 



±s(t) = T](t), 



(13) 



and the probability distribution of the endpoints of paths starting at a certain initial point is given by 
the path integral 



P(x b t b \x a ta) 



Drj / T>x exp 



dtH(rj(t)) 



5[x - 



(14) 



with the initial condition Xs(t a ) — x a . The final point is, of course, Xb — xs{tb)- The function H(rj) is 
defined by the negative logarithm of the non-Gaussian distributions D(x), such that 



e -H{x) = jj^ 



(15) 



The path integral is defined a la Feynman [j35f| by slicing the time axis at times t n = en with n = 0, . . . , N, 
and integrating over all x(t n ). At the end, the limit N is taken. In this way, we select from the space 
of all fluctuating paths a well-behaved set of measure zero, which is in most cases sufficient to obtain 
the correct limit N — > oo, in particular it will be for the typical non-Gaussian distributions encountered 
in stock markets. This is completely analogous to Riemann's procedure of defining ordinary integrals, 
which are approximated by sums over values of a function on a set of measure zero. 

With the help of the path integral (fTi]), the correlation functions of the noise in the path integral 
( |l4| ) can be found by straightforward functional differentiation. For this purpose, we express the noise 

distribution P[rj\ = exp — J^ b dt H(rj(t)) as a Fourier path integral 



/'['/] = / P'll f^exp 



dtMt)n(t)-H(p(t))} 



and note that the correlation functions can be obtained from the functional derivatives 



Vp 
2^ 



Sp(h) 5p(t n ) 



S i f tb dtp(t)r)(t) 

e Jt a 



f t b dtH{p{t)) 



After n partial integrations, this becomes 

- f t b dtH(p(t)) 



5 - f tb dtH(p(t)) 

e Jt " 



5p{t n ) 



Sp(ti) Sp(t n ) 



p(t)=0 



- f b dtH(p(t)) . 



(16) 



(17) 



(18) 



By expanding the exponential e 



in a power series using ( ^0| ) , it is straightforward to calculate 



dtH(rj(t)) 



Cl, 



WhMh)) = Z- 1 I Vrjriih)^) exp 
= c 2 S(h - t 2 ) + cl, 
(V(h)v(t2)v(h)) = Z- 1 [ P7 ? ry(i 1 )7 ? (t 2 )7 7 (t 3 )exp 



dtH(r)(t)) 



(19) 



(20) 



dtH(rj(t)) 
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= caS(t! - h)8(ti - t 3 ) 

+ c 2 ci [<y(i! - t 2 ) + 5{t 2 - i 3 ) + <y(ti - i 8 )] , 

(v{h)r}{t2)r}{h)v{U)) = Z- 1 [ V v v(tiMhMh)v(U) exp 



(21) 



dtH(r)(t)) 

= Ci5{h - t 2 )5{h - t 3 )S(h - U) 

+ C3C1 [S(ti — t 2 )5(ti — i 3 ) + 3 cyclic perms] 

+ c 2 <?\ [5{ti — t%) + 5 pair terms] 

+ Ac\ [5{h - t 2 )S(t 3 - i 4 )<K*i - h)5(t 2 - t A )S{t! - t 4 )S(t 2 - t 3 )] . 



(22) 



where Z is the normalization integral 

Z = I T)r\ exp 



dtH(r)(t)) 



(23) 



The higher correlation functions are obvious generalizations of these equations. The different contribu- 
tions on the right-hand side of Eqs. (p^|)-( |22|) are distinguishable by their connectedness structure. 
An important property of the probability (|14|) is that it satisfies the semigroup property of path integrals 



P(x c t c \x a t a ) = / dxb P(x c t c \x b t b )P(x b t b \x a t a ). 



(24) 



The the experimental asset distributions do satisfy this property approximately For truncated Levy 
distributions this is shown in Fig. [j]. 
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FIG. 1. Left-hand side: Cumulative distributions obtained by repeated convolutions of the 15-min distribution 
(from Ref. |B|). Apart from the tails, the semigroup property (^) is reasonably well satisfied. Right-hand side: 
Falloff of kurtosis is slower than — t a ) expected from convolution property (see Ref. 

Apart from the far ends of the tails, the semigroup property (G4) is reasonably well satisfied. 



B. Levy-Khintchine Formula 

If is sometimes useful to represent the Hamiltonian in the form of a Fourier integral 

Hip) = - J dxe ipx F(x). (25) 
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Due to the the special significance of the linear term in H (p) governing the drift, this is usually subtracted 
out of the integral by rewriting ( ^5|) as 

ftM-*, + /*(*--i-*.) W M 

The first subtraction ensures the property H r (0) = which guarantees the unit normalization of the 
distribution. This subtracted representation is known as the Levy-Khintchine formula, and the function 
F(x) is the so-called of the distribution. Some people also subtract out the quadratic term and write 

2 p 

H r (p) = irp + yp 2 + dx (e ipx - 1) F(x). (27) 

and use a function F(x) which has no first and second moment ( i.e., J F(x)x — 0, J F(x)x 2 — 0. to 
avoid redundancy in the representation. 



C. Fokker-Planck-Type Equation 



The J-functional may be represented by a Fourier integral leading to 



P{x b t b \x a t a ) = / Vrj 



Vp 
2^ 



exp 



dt 



ip(t)x s (t) - ip(t) v (t) - H( v (t)) 



Integrating out the noise variable rj (£) at each time t we obtain 



P(x b t b \x a t a ) 



Vp 
2^ 



exp 



dt[ip(t)x s {t)-H(p(t))} 



(28) 



(29) 



Integrating over all xs(t) with fixed end points enforces a constant momentum along the path, and we 
remain with a single integral 



dp 

P{x b t b \x a t a ) = I — exp [-(4 - t a )H{p) + ip(x b - x a )\ 

Z7T 



(30) 



From this representation it is easy to verify that this probability satisfies a Fokker-Planck-type equation 

d t P(x b t b \x a t a ) = -H(-id x )P{x b t b \x a t a ). (31) 

The general solution ip(x, i)of this differential equation with the initial condition tf)(x, 0) is given by the 
path integral generalizing (|1J 



tp{x, t) — T>rj exp 



dtH(r)(t)) 



iplx- dt'rj{t') 



(32) 



To verify that this satisfies indeed the Fokker-Planck-type equation ( |3l| ) we consider tp(x,t) at a slightly 
later time t + e and expand 



ip(x, t + e) 



T>7] exp 
Vrj exp 



dtH(r)(t)) 



dtH{r/{t)) 



ip yx — 
ip I x 



dt'r)(t') - 
t 

d£j]{t' 



t+e 



dt'r)(t') 



i>' (x - J dt'ri(t'fj £ 6 dt'r](t') 



3!^ ' X 



dt'r/{t' 



dt'r)(t' 



t+e 



dt 1 dt 2 ri(t 1 )T){t2) 



(33) 



L+e 



dtidt 2 dt 3 T7(ti)77(i 2 )/7(i3) 



+ ^V (4) ( x ~ J J dhdhdtsdU vihMhMhWh) 
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Using the correlation functions ([L9l)— (E2I) we obtain 



ifj(x, t + e) 



T>r\ exp 



dtH(r)(t)) 



-ec x d x + (ec 2 + e 2 Cl ) i<9 2 - (ec 3 + 3e 2 c 2 c 1 ) -8% 



(34) 



+ (ec 4 + e 2 4c 3 ci + e 2 3c 2 + e 3 c 2 c 2 + e 4 c?)^ 



ip x 



dt'r](t') 



In the limit e — > 0, only the linear terms in e contribute, which are all due to the connected parts of 
the correlation functions of rj(t). The differential operators in the brackets can now be pulled out of the 
integral and we find the differential equation 



d t ip(x,t) 



-cid x + c 2 i<9 2 - c 3 -d* + c A j { di 



ip(x,t) . 



(35) 



We now replace c\ — > r x and express using ( |To| ) the differential operators in brackets as Hamiltonian 
operator —H rsc (—id x ). This leads to the Schrodinger-like equation 



d t ip(x,t) 



-id x )ip(x,t) . 



(36) 



Due to the many derivatives in H(id x ), this equation is in general non-local. This can be made explicit 
with the help of the Levy-Khintchine weight function F(x) in the Fourier representation (|2^). In this 
case, the right-hand side 



H(id x )ijj(x,t) = / dx' e- x ' d * F(x')ip(x,t) = dx' F(x')ijj(x - x' ,t) 



(37) 



and the Fokker-Planck-like equation ( |36| ) takes the form of an integral equation. Some people like to use 
the subtracted form ( ^7j ) of the Levy-Khintchine and arrive at the integro-differential equation 



d t tp(x,t) 



-c\d x 



?9 2 
2 



ip(x,t)+ / dx' F(x')ip(x - x',t). 



(38) 



The integral term can then be treated as a perturbation to an ordinary Fokker-Planck equation. 

By a similar procedure as in the derivation of Eq. (35) it is possible to derive a generalization of Ito's 
rule (^|) to functions of noise variable with non-Gaussian distributions. As in (^) we expand f{x{t + e)): 



f(x(t + e)) = f(x(t))+ f'(x(t)) / dt'x(t') 



t+e 



;/"(*(*)) 



i+<- 



dti 



t+e 



dt 2 x(ti)i(t 2 



(39) 



+ ^/ (3) (^W) / dt x 



t + c 



dU 



t+e 



t 3 x(ti)x(t 2 )x(t 3 ) + 



where x{t) = r](t) is the stochastic differential equation with a nonzero expectation value (r](t)) = c\. In 
contrast to (|^) which had to be carried out only up to second order in x, we must now keep all orders 
in the noise variable. Evaluating the noise averages of the multiple integrals on the right-hand side using 
the correlation functions (|l9|)-(f22"|), we find the time derivative of the expectation value of an arbitrary 
function of the fluctuating variable x(t) 



(f(x(t + e))) = (f(x(t))) + (f(x(t)))e Cl + -(/"(x(t)))(ec 2 +e 2 c?) 

+ ^</ (3) (^))>(ec3+ e 2 c 2 c 1+ e 3 cl) + . . 



'Cid x + c 2 ^d 2 x - c^dl 



(f(x(t)))+0(e 2 ). 



(40) 
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After the replacement c\ — ► r x the function f(x(t)) obeys therefore the following equation: 

(f(x(t)))=-H rx (id x )(f(x(t))). (41) 

Taking out the lowest-derivative term this takes a form 

(f(x(t))) = (d x f(x(t)))(x(t)) - H rx (id x )(f(x(t))). (42) 

In postpoint time slicing, this may be viewed as the expectation value of the stochastic differential 
equation 

/(*(*)) = d x f(x(t))x(t) - H rx (id x )f(x(t)). (43) 

This is the direct generalization of Ito's rule (||). 

For an exponential function f(x) = e Px , this becomes 

—e Px{t) = \x(t) - H r (iP)} e Px W. (44) 
at 

As a consequence of this equation for P = 1, the rate rs with which a stock price S(t) = e x ^ grows 
according to formula (|^) is now related to r xs by 

rs = r xs - H(i) = r xs - [H(i) - iH'(0)] - -H r „ g (i), (45) 

which replaces the simple Ito relation rs — r xs + a 2 /2 in Eq. Recall the definition H(p) = H(p) — 
H'(Q)p in Eq. (12). The corresponding generalization of the left-hand part of Eq. (0) reads 

| = x(t) - H(i) = ±{t) - [H{i) - iH'(0)} = ±{t) - r xs - H Txs (t). (46) 
The forward price of a stock must therefore be calculated with the generalization of formula (^): 

(S(t)) = 5(0)e rst = S(0)(e r * st+ fo dt ' = S(0)e~ H ^s^ = S (0)e^s-[H(i)-iH'(o)]t} ^ (47) 

Note that may derive the differential equation of an arbitrary function f(x(t)) in Eq. ( f4l| ) from a simple 
mnemonic rule, expanding sloppily 

f(x(t + dt)) = f(x(t) + xdt) = f(x(t)) + f'(x(t))x(t)dt + ^}"{x{t))x 2 (t)dt 2 

+ ^f {3) ^m 3 (t)dt 3 + ... , (48) 

and replacing (^]), 

(x(t))dt -> cidt, (x 2 {t))dt 2 -> cadi, (x 3 {t))dt 2 -> c 3 di, . . . . (49) 



IV. MARTINGALES 

In financial mathematics, an often-encountered concept is that of a martingale J37j. The name stems 
from a casino strategy in which a gambler doubles his stake each time a bet is lost. A stochastic variable 
is called a martingale, if its expectation value is time-independent. The noise variable rj(t) with vanishing 
average is a trivial martingale. 
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A. Gaussian Martingale 

For a harmonic noise variable, the exponential e^o dt n ^ t ' a is a nontrivial martingale, due to Eq. (^). 
For the same reason, a stock price S(t) = e x ^ with x(t) obeying the stochastic differential equation ([l ) can 
be made a martingale by a time-dependent multiplicative factor the quantity e~ rst e x = e~ r ^s t e x ~ a t / 2 . 
An explicit distribution which makes Sit) = e x ^ a martingale is 

}M/ i |„ j. -\ e / i X b 2-a r xs(^b ^a)] 



P M {x b t b \x a t a ) = - =exp - L " " , sv \ a/J ^ (50) 

I can easily be verified by direct integration that the expectation value of S(t) = e xtyt ^ is time-independent 
with this distribution: 

(S(t b )) = (e x ^) = [ dx b e Xb P M (x b t b \x a t a ) (51) 



is independent of the time t b . 

At this place we can make an important observation: There exists an entire family of distributions for 
which S(t) — e x ^ is a martingale, namely 

P Mr (x b t b \x a t a )^ - e ~ rt ™ P {- [Xb - X "-J x{tb - ta)]2 \. (52) 

for any r and r x = r — a 2 /2. Such distributions which differ by the drift r are called equivalent. The 
prefactor e~ rt is referred to as a discount factor with the rate r. 

B. Non-Gaussian Martingales 

For S(t) = e x &' with an arbitrary non-Gaussian noise r)(t), there are many ways of constructing 
martingales, the relation ( ]45| ) allows us to construct immediately the simplest martingale 

1. Natural Martingale 

An obvious generalization of the Gaussian expression e-^o dt vt - 1 -* a t ^ 2 is the exponential 
e I dt )+H(i)t ^ w j 10se expectation value is time-independent due to Eq. ([fo]). Expressed in terms of 
the associated stock price S(t) = e xtyt \ we obtain the simplest martingale 

e~ rst S{t) = e - rst e r *s t+ fo dt '" ( *' ) (53) 

if rs and r Xs are related by rs = r Xs — H(i) — —H Xs (i). 

It is easy to write down a distribution function which makes S(t) — e x<yt ^ itself a martingale: 

P M (x b t b \x a t a ) = e~ rst [vn ( Pxcxpj- / " dtH r Mt))\ 6[x - rj\. (54) 



There exists also here an entire family of equivalent distribution functions for which the integrals of the 
type ( |5l| ) over e Xb are independent of t b . This is done as follows. One obvious family of this type is a 
straightforward generalization of the Gaussian family (|52|) which we shall refer to as natural martingales: 



P Mr (x b t b \x a t a )=e- rt / V V / Vxexpl - / dt H r .(»/(t)) \ S[x - rj\, (55) 



with arbitrary r = r x — H{i) = —H,, x (i). Indeed, multiplying this with e Xh and integrating over x b gives 
rise to a 5- function S(p — i) and produces the same result e Xa for all times t b . The path integral is solved 
by the Fourier integral [compare (|3C|)1 

P Mr (x b t b \x a t a ) = e~ rt [ ^ exp [ip(x b - x a ) - (t b - t a )H rx (p)\ • (56) 
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2. Esscher Martingale 



In the literature on mathematical finance, much attention is given another family of equivalent mar- 
tingale measures. It has been used a long time ago to estimate risks of actuaries | p8 fl and introduced 
more recently into the theory of option prices [p9|j4C| 1 where it is now of wide use 0P|||. This family 
is constructed as follows. Let D(x) be an arbitrary distribution function with a Fourier transform 



D(x) 



2ir 



(57) 



and -ff(O) = 0, to guarantee a unit normalization J dxD(x) — 1. We now introduce an Esscher- 
transformed distribution function. It is obtained by slightly tilting the initial distribution D{x) by 
multiplication with an asymmetric exponential factor e 6x : 



D e (x) = e H ^ e 0x D(x). 



(58) 



The constant prefactor e H ^ l6%> is necessary to conserve the total probability. This distribution can be 
written as a Fourier transform 



2tt 



D\x) -- 

with the Esscher-transformed Hamiltonian 

H e (p) = H(p + i8) - H{i 



(59) 



(60) 



Since H 6 (0) = 0, the transformed distribution is properly normalized J dxD e (x) — 1. We now define the 
Esscher-transformed expectation value 



(F{x)) b 



It is related to the original expectation value by 



dxD e (x)F{x). 



(F(x)) 6 = e H{ie) (e 8x F(x)) 
For the specific function F(x) — e x , Eq. ( |6^ ) becomes 

(e x ) e = er H<>{l) = e H " (iB) (e {e+1)x ) = e H 9 (i0)-H 9 (i0+i) _ 



(61) 



(62) 



(63) 



Applying the transformation ((59|) to each time slice in the general path integral (|TJ), we obtain the 
Esscher-transformed path integral 



dt 



6r,(t) - H r ( V (t)) \6[±-ri\, 



(64) 



The solution is given by the Fourier integral [compare ( p6| 



P e (x b t b \x a t a ) = e- rst e H ^s ^ I ^ exp ip(x b - x a ) - (t b - t a )H r (p + i6 

-,2.711 ° 



Let us denote the expectation values calculated with this probability by 
S{t) = e x ^ the time dependence 



(Sit))' 



(65) 

Then we find for 
(66) 



This equation shows, that the exponential of a stochastic variable x(t) can be made a martingale with 
respect to any Esscher-transformed distribution if we remove the exponentially growing factor exp(r e t) 
with 
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r° ee -H° mg (<) = -H r9B (i + iff) + ff, xs (i6). (67) 

Thus, a family of equivalent martingale distributions for the stock price S(t) — e x<yt ^ is 

P Me {x b t b \x a t a ) EE e- rst P e (x b t b \x a t a ) (68) 

for any choice of the parameter 9. 

For a harmonic distribution function (p"o|), the Esscher martingales and the previous ones are equivalent. 
Indeed, starting from (5C|) in which r$ = r xs + cr 2 /2, the Esscher transform leads us after a quadratic 
completion to the family of natural martingales (pOI) with the rate parameter r — r xs + da 2 . 

3. Other Non-Gaussian Martingales 

Many other non-gaussian martingales have been discussed in the literature. Mathematicians have 
invented various sophisticated criteria under which one would be preferable over the others for calculating 
financial risks. Davis has introduced a so-called utility function fl33| which is supposed to select optimal 
martingales for different purposes. 

For the upcoming development of a theory of option pricing, only the initial natural martingale will 
turn out to be relevant. 



V. OPTION PRICING 

The most important use of path integrals in financial markets is made in the determination of a 
fair price of financial derivatives, in particular options. Options are an ancient financial tool. They 
are used for speculative purposes or for hedging major market transactions against unexpected changes 
in the market environment. These can sometimes produce dramatic price explosions or erosions, and 
options are supposed to prevent the destruction of huge amounts of capital. Ancient Romans, Grecians, 
and Phoenicians traded options against outgoing cargos from their local seaports. In financial markets, 
options are contracted between two parties in which one party has the right but not the obligation to do 
something, usually to buy or sell some underlying asset. Having rights without obligations has a financial 
value, so option holders must pay a price for acquiring them. The price depends on the value of the 
associated asset, which is why they are also called derivative assets or briefly derivatives . Call options 
are contracts giving the option holder the right to buy something, while put options entitle the holder 
to sell something. The price of an option is called premium. Usually, options are associated with stock, 
bonds, or commodities like oil, metals or other raw materials. In the sequel we shall consider call options 
on stocks, to be specific. 

Modern option pricing techniques have their roots in early work by Charles Castelli who published in 
1877 a book entitled The Theory of Options in Stocks and Shares. This book presented an introduction 
to the hedging and speculation aspects of options. However, it still lacked a sound theoretical base. 
Twenty three years later, Louis Bachelier offered the earliest known analytical valuation for options in 
his dissertation at the Sorbonne . It is curious that Bachelier discovered the treatment of stochastic 
phenomena five years before Einstein's related but much more famous work on Brownian motion p8| , and 
23 years before Wiener's mathematical development [Q. The stochastic differential equations considered 
by him still had an important defect of allowing for negative security prices, and for option prices exceeding 
the price of the underlying asset. Bachelier's work was continued by Paul Samuelson, who wrote in 1955 
an unpublished paper entitled Brownian Motion in the Stock Market. During that same year, Richard 
Kruizenga, one of Samuelson's students, cited Bachelier's work in his dissertation Put and Call Options: 
A Theoretical and Market Analysis. In 1962, another dissertation, this time by A. James Boness, focused 
on options. In his work, entitled A Theory and Measurement of Stock Option Value, Boness developed a 
pricing model that made a significant theoretical jump from that of his predecessors. More significantly, 
his work served as a precursor to that of Fischer Black and Myron Scholes, who in 1973 introduced their 
famous Black and Scholes Model IpCf which, together with the improvements introduced by R. Merton, 
earned the Nobel prize in 1997.Q 



1 For F. Black the prize came too late — he had died two years earlier. 
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A. Black-Scholes Option Pricing Model 



In the early seventies, Fisher Black was working on a valuation model for stock warrants and observed 
that his formulas resembled very much the well-known equations for heat transfer. Soon after this, Myron 
Scholes joined Black and together they discovered an approximate option pricing model which is still of 
wide use. It is an improved version of a previous model developed by A. James Boness in his Ph.D. thesis 
the University of Chicago. 

The Black and Scholes Model is based on the following assumptions: 

1. The logarithm of the returns is normally distributed. We remarked before that there are considerable 
deviations which call for improvement of the model to be developed below. 

2. Markets are efficient. This assumption implies that the market operates continuously with share 
prices following a continuous stochastic process without memory. It also implies that different 
markets have the same asset prices. 

This is not quite true. Different markets do in general have slightly different prices. Their differences 
are kept small by the existence of arbitrage dealers. There also exist correlations over a short time 
scale which make it possible, in principle, to profit without risk from statistical arbitrage. This 
possibility is, however, strongly limited by transaction fees. 

3. No commissions are charged. 

This is not true. Usually market participants have to pay a commission to buy or sell assets. Even 
floor traders pay some kind of fee, although this is usually very small. The fees payed by individual 
investors is more substantial and can distort the output of the model. 

4. Interest rates remain constant and known. 

The Black and Scholes model assumes the existence of a risk-free rate to represent this constant 
and known rate. In reality there is no such thing as the risk-free rate. As an approximation, one 
uses the discount rate on U.S. Government Treasury Bills with 30 days left until maturity. During 
periods of rapidly changing interest rates, these 30 day rates are often subject to change, thereby 
violating one of the assumptions of the model. 

5. The stock pays no dividends during the option's life. dividends 

Most companies pay dividends to their share holders, so this is a limitation to the model since 
higher dividend lead to lower call premiums. There is, however, a simple possibility of adjusting 
the model to the real situation by subtracting the discounted value of a future dividend from the 
stock price. 

6. European exercise terms are used. European exercise terms imply the exercise of an option only 
on the expiration date. This is in contrast to the American exercise terms which allow for this 
at any time during the life of the option. This greater flexibility makes an American option more 
valuable than the European one. 

The difference is, however, not dramatic in praxis because very few calls are ever exercised before 
the last few days of their life, since an early exercise means giving away the remaining time value 
on the call. Different exercise times towards the end of the life of a call are irrelevant since the 
remaining time value is very small and the intrinsic value has a small time dependence, barring a 
dramatic event right before expiration date. 

Since 1973, the original Black and Scholes Option Pricing Model has been the subject of much attention. 
In the same year, Robert Merton included the effect of dividends. Three years later, Jonathan 
Ingerson relaxed the assumption of no taxes or transaction costs, and Merton removed the restriction 
of constant interest rates. At present we are in a position of being able to determine quite reliably the 
values of stock options. 

The relevance of path integrals to this field was recognized first in 1988 by a theoretical physicist J.W. 
Dash, who wrote two unpublished papers on the subject entitled Path Integrals and Options I and 77 
p2^| . Since then many theoretical physicists have entered the field, and papers on this subject have begun 
appearing on the Los Alamos server . 
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B. Evolution Equations of Portfolios with Options 



The option price 0{t) has a larger fluctuations than the associated stock price. It usually varies with 
an efficiency factor dO(S(t),t)/dS(t). For this reason it is possible, in the ideal case of Gaussian price 
fluctuations, to guarantee a steady growth of a portfolio by mixing Ns(t) stocks with No{t) options and 
a certain amount of cash which is usually kept in the form of bonds, whose number is denoted by N B (t). 
The composition [ATg(f), No(t), N B (t)] is referred to as the strategy of the portfolio manager. The total 
wealth has the value 



W{t) = N s {t)S{t) + N o (t)0{S,t) + N B (t)B{t). (69) 

The goal is to make it grow with a smooth exponential curve without fluctuations 

W(t) w r w W(t). (70) 

As we shall see immediately, this is possible provided the short-term bonds (usually those with 30 days 
to maturity) grow without fluctuations. 

Bit) « r B B{t). (71) 

Under this assumption, the rate r B is referred to as riskfree interest rate This assumption is fulfilled in 
true markets only approximately since there can always be events which change the value of short-term 
bonds excessively. 

The existence of arbitrage dealers will ensure that the growth rate rw is equal to that of the short-term 
bonds 



rw ~ tb- (72) 



Otherwise the dealers would change from one investment to the other. 
In the decomposition (169), the desired growth (JTOj) reads 



N s {t)S(t) + N o (t)0(S,t) + N B (t)B(t) + N s (t)S(t) + N o (t)0(S,t) + N B (t)B{t) 

= r w [N s (t)S(t) + N o (t)0(S, t) + N B (t)B(t)} . (73) 



Due to ( ]7l| ) and (|72[), the terms containing N B (t) without a dot drop out. Moreover, if no extra money 
is inserted into or taken from the system, i.e., if stocks, options, and bonds are only traded against 
each other, this does not change the total wealth, assuming the absence of commissions. This so-called 
self-financing strategy is expressed in the equation 

N s (t)S(t) + N (t)0(S,t) + N B (t)B(t) = 0. (74) 



Thus the growth equation (70) translates into 

W(t) = N S S + N + N B B = r w (n s S + N + N B B^j . (75) 

Due to the equality of the rates rw = r B and Eq. (|7l|), the entire contribution of B(t) cancels, and we 
obtain 

N S S + N o = r w (N S S + N 0) (76) 

The important observation is now that there exists an optimal ratio between the number of stocks N$ 
and the number of options No, which is inversely equal to the efficiency factor 

N 3 (t) _ dO(S(t),t) 

N (t) dS(t) ■ 1 ' 

Then Eq. (|76| ) becomes 

N s S + N d = r w f-^ + o) N . (78) 
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The two terms on the left-hand side are treated as follows: First we use the relation ( |77j ) to rewrite 

dO(S,t) 



N S S 



dS 



-S 



and further, with the help of Eq. ( |46[) , as 



N S S = 



dO{S, t) 
dx 



dO{S,t) S 
dx S' 



[x - H{i) 



(79) 



(80) 



In the second term on the left-hand side of (|78|) , we expand the total time dependence of the option price 
in a Taylor series 



dO 
~dt 



0(x(t) + x{t) dt, t + dt)- 0(x{t), t) 



1 

~dt 

80 dO 
dt dx 



1 d 2 Q 

2 9^" 



: 2 dt 



1 d 3 
3! a^" 



; 3 dt 2 



(81) 



We we have gone over to the logarithmic stock price variable x(t) rather than S{t) itself. Some of 
the derivatives on the right-hand side are denoted by special symbols in financial mathematics: the 
quantities = dO/dt, A = dO/dS = dO/dxS, and T = d 2 0/dS 2 = (d 2 0/dx 2 - dO/dx)/S 2 are called 
the "Theta" , "Delta" , and "Gamma" of the option. Another derivative with a standard name is the 
"Vega" V = dO /da. 

The expansion ( |Sl| ) is carried to arbitrary powers of x as in (|48|). It is, of course, only an abbreviated 
notation for the proper expansion in powers of a stochastic variable to be performed as in Eq. (|39|). 
Inserting (pTI) and (18^) on the left-hand side of Eq. (|78|) , this becomes 



N S S + N 



No^- [±-H[i)} 



+ N 

= N 



dx 

dO 
~dt 



dO 

dx 



-...SO do 



ld 2 0. 2l ISO.o, 

i9o. 2l iao.,, 



(82) 



Remarkably, the fluctuating variable x drops out in this equation, which therefore looses its stochastic 
character. At the same time, it becomes independent of the growth rate rg of the stock price. This is 
the reason why the total wealth W(t) increases without fluctuations. 
We now treat the Taylor series 



ld 2 



2 dx 2 



x 2 dt + -^ri 3 dt 2 



6 dx 3 



in the same way as the expansion (48), using the rules 



ld 2 



■ 2 dt 



ld 3 



: 3 dt 2 - 



|), such that (|83j) becomes 
= -H(id x )0. 



2 dx 2 6 dx 3 

In this way we find for the option price 0{x,t) the Fokker-Planck-like differential equation 

dO dO 

— = r w O - r xw — + H{id x )0, 

where we have defined, by analogy with ([i"5|), an auxiliary rate parameter 

r xw = r w + H(i). 



(83) 



(84) 



(85) 



(86) 



Note, however, that in contrast to the relation between r xs and r$ defined for a fluctuating stock price 
S(t) and its logarithm x(t), the parameter r xw does not have the physical interpretation of governing 
the logarithm of W(t) since the absence of fluctuations in the wealth W(t) makes log W(t) grow linearly 
with the same rate rw that governs the exponential growth e rwt of W{t) itself. 
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If we rename t as t a , The general solution of the differential equation (pq), which at some time t — tb 
starts out like 8(x — x b ), has the Fourier representation 

P(x b t b \x a t a ) = e -^(*»-*«) y°° !pp(*>-*0 exp {- [H(p) + ir Slv p] (t b - t a )} , (87) 

if the initial variables x a and t a are identified with x and i, respectively. A convergent integral exists 
only for t b > t a . 

Comparing this expression with Eq. fl56|) we recognize it as a member of the family of equivalent 
martingale measures P Mr (xbtb\x a t a ), in which the discount factor r coincides with the risk-free interest 
rate rw 



C. Option Pricing for Gaussian Fluctuations 

For Gaussian fluctuations where H(p) = a 2 p 2 /2, the integral in (|87| ) can easily be performed and yields 

P{x b tb\x a t a ) = 6(4 - t a ) e~ rw(tb - ta) exp / - ^ - (*6 - *a)]' 



This probability distribution is obviously the solution of the path integral 



P(x b tb\x a t a ) = e(tb-t a )e- rw ^- t ^ j Pxexpj-^ / " [x - r Xw ] 2 \ . (89) 



Recalling the discussion in Section IV, the distribution function ( |88| ) is recognized as a member of the 
equivalent family of martingale distributions J52|) for the stock price S(t) = e x ^\ It is the particular 
distribution in which the discount factor contains the risk-free interest rate rw, i-e., ( p8[ ) is equal to the 
martingale distribution P Mrw (xbtb\x a t a )- This distribution is referred to as the risk-neutral equivalent 
martingale distribution. 

An option is written for a certain strike price E of the stock. The value of the option at its expiration 
date t b is given by the difference between the stock price on expiration date and the strike price: 

0(x b ,t b ) = 0(S b - E)(S b -E) = 0(x b - x E ){e x » - e XE ), (90) 

where 

x E = \ogE. (91) 

The Heavisidc function accounts for the fact that only for Sb > E it is worthwhile to execute the option. 
From (Q) we calculate the option price at an arbitrary earlier time using the time evolution amplitude 

(ID 

/CO 
dx b 0(x b ,t b ) P Mrn '(x b t b \x a t a ). (92) 
-oo 

Inserting (Q) we obtain the sum of two terms 

0(x a ,t a ) = Os(x a ,t a ) - E (x a , t a ), (93) 

where 

n, *\ e-^to-V f°° , f [xb-x a -r xw (tb-t a )] 2 \ 

O s {x a ,t a ) = = / dx b expixb — >. 94 

y/2Tr<T 2 (t b -t a ) J XE I 2a\t b -t a ) J 

and 

\x b X a Txw {tb t a )] 



E (x a ,t a ) ^Ee- rw ^- 1 ^— = } t f dx b exp{- , 

y/2ira*(jk b -t a ) J XE I 2a 2 (t b -t a ) 



(95) 
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In the second integral we set 



X- = x a + r Xw (t b - t a ) = x a + ( r w - -cr ) (t b - t a ), (96) 



and obtain 



E (x ai t a ) = E dx b exp\- - b . (97) 

y/2ncr 2 (t b - t a ) J XE - X _ { 2a 2 {t b -t a )) 

After rescaling the integration variable x b — > — £<7a/^6 — t a , this can be rewritten as 

E (x a ,t a ) = e-^^-^EN^), (98) 

where N(y) is the Gaussian distribution function 

N (y)= f V JL e -f/\ (99) 



2tt 

evaluated at 

_ x--x E \og[S(t a )/E]+r xw (t b -t a ) 



y- 



V° 2 (ta-t b ) V° 2 (ta-t b ) 

\og[S(t a )/E] + (rw - \° 2 ) (t b - t a ) 
^e 2 {t a -t b ) 



(100) 



The integral in the first contribution ( |94| ) to the option price is found after completing the exponent 
in the integrand quadratically as follows: 



_ [x b - x a - r Xw (t b - t a )} 2 
X " 2a 2 (t b -t a ) 

_ [x b -x a - (r Xw + a 2 )(t b - t a )} 2 - 2r w a 2 (t b - t a ) - 2x a a 2 (t b - t a ) 

2a 2 (t b -t a ) 

Introducing now 

x + =x a + (r Xw + a 2 ) (t b - t a ) = x a + yr w + ^cr 2 ^ (t b - t a ), 



(101) 



(102) 



and rescaling x b as before, we obtain 

O s (x a ,t a ) =S(t a )N(y+), (103) 

with 

x + -xe _ \og[S(t a )/E] + (r xw + a 2 ) (t b - t a ) 



y+ 



V° 2 (ta-t b ) V° 2 (ta-t b ) 

log[S(t a )/E] + (r w + \a 2 ) (t b - t a ) 



(104) 



y/<J 2 (t a -t b ) 
The combined result 

0(x a ,t a ) = S(t a )N(y + ) - e- rw( - tb -^E N(y_) (105) 

is the celebrated Black-Scholes formula of option pricing. 

In Fig. H we illustrate how the dependence of the call price on the stock price varies with different 
times to expiration t b — t a and with different volatilities a. 
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20 40 60 



FIG. 2. Left: Dependence of call price O on the stock price S for different times before expiration date (in- 
creasing dash length: 1, 2, 3, 4, 5 months), from left to right: 1, 2, 3, 5, 6 months). The parameters are E — 50 US$, 
a = 40%, rw = 6% per month. Right: Dependence on the strike price E for fixed stock price 35US$ and the 
same times to expiration (increasing with dash length). Bottom: Dependence on the volatilities (from left to 
right: 80%, 60%, 20%, 10%, 1%) at a fixed time tt — t a = 3 months before expiration. 



D. Option Pricing for Non-Gaussian Fluctuations 



For non-Gaussian fluctuations, the option price must be calculated numerically from Eqs. (|92j) and 
(pp|). Inserting the Fourier representation (87) and using the Hamiltonian 



H r xvv (p) = Hip) + ir XwP 



(106) 



defined as in (ff3), this becomes 



Oix a ,t a ) = / dx b (e x > - e XE )P{x b t b \x a t a ) 



_ p -r w {t b -t a ) 



dxb (e 



dP e ip(x b -x a )-H rxw (p)(t b -t a ) 

2tt 



(107) 



The integrand can be rearranged as follows: 



£ X a e (ip+l)(x b -X a ) _ p XE p ip(x b -Xa) 



(108) 



Two integrations are required. This would make a numerical calculation quite time consuming, Fortu- 
nately, one integration can be done analytically. For this purpose we write the integral in the form 



0(x a ,t a ) = e-^to-**) 



dxh 



QP r ip(Xh-Xg) 

2tt 



' " aj f(x a ,x E ;p), 



(109) 



with 



(110) 



We have suppressed the arguments x a , xe, % ~ ta m fip), for brevity. The integral over x b in (109) runs 
over the Fourier transform 



fix b -x a ) = 



dP r ip(x b -x a ) 

2tt 



f(p), 



(111) 
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of the function f(p). It is then convenient to express the integral J^dxb in terms of the Heaviside 
function Q(x b — xe) as dxb ®{xb — xe) and use the Fourier representation 

Q(x b -x E )= I p.^—e-^-^). (112) 
of the Heaviside function to write 

' dx b f(x b -x a ) = r dx b [ p—L^e-W't-'^fixb-Xa). (113) 



Inserting here the Fourier representation (111), we can perform the integral over x b and obtain the 
momentum space representation of the option price 

°^ u)= Lt^' n) jhi m (114) 

For numerical integrations, the singularity at p = is inconvenient. We therefore use the well-known 
decomposition 

= -+nS(p), (115) 



p + 17] p 

to write 

2 7_oo 2tt p 

We have used the fact that the principal value of the integral over 1 jp vanishes to subtract the constant 
/(0) from e l P( x E-x a ) J(p). After this the integrand is regular, does not need any more the principal-value 
specification, and allows for a numerical integration. 

For x a very much different from xe, we may approximate 

■"•''*'-'- >f{ » } - m ~U*«-*zW». (117) 



, 2tt p 2 

where e(x) = 1 + 20(x) is the step function, and obtain 

0(x a ,t a ) wi[l + e(x a - x E )] /(0). (118) 

Using ( |86| ) we have e~ Hrx w W _ e rw(t b -t a ) ^ an( j gmce e ~ Hr tc W (°) _ we see ^h a t (9(a; Qj i a ) goes to zero 
for x a — > — oo and has the large- ir a behavior 

0(x a ,t a ) » (e x ° - e ^ e -'-»'(t i ,-ia)^ = _ e -rw(t„ ( 119) 



This is the same behavior as in the Black-Scholes formula ( |105| ). 

In Fig. ^ we display the difference between the option prices emerging from our formula ( |116| ), with a 
truncated Levy distribution of kurtosis k = 4, and the Black-Scholes formula ( |105| ) for the same data as 
in the upper left of Fig. 0. 



0.4 0(S)-0 BS (S) 

0-2 // / / V\C\ 



20 40, 



-60 



-0.2 
-0.4 



FIG. 3. Difference of the call price O obtained from truncated Levy distribution with kurtosis n — 4 and the 
Black-Scholes price as a function of the stock price S for different times before expiration date (increasing dash 
length: 1, 2, 3, 4, 5 months). The parameters are E = 50 US$, a = 40%, rw = 6% per month. 
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VI. CONCLUSION 



The stochastic calculus and the option pricing formulas developed in this paper will be useful for 
estimating financial risks of a variety of investments. In particular, it will help developing a more realistic 
theory of fair option prices. More details can be found in the textbook Ref. p5[. 
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